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FOREWORD 


This  interim  report  summarizes  the  experimental  results  on  OH  LIF 
Temperature  Measurements  to  High  Pressure  Solid  Propellant  Flames  in  the 
Combustion  Research  Laboratory  at  the  Air  Force  Astronautics  Laboratory 
(AFAL),  Edwards  Air  Force  Base  CA.  The  report  period  is  from  January  1988  to 
January  1989.  The  AFAL  project  manager  was  Tim  Edwards. 
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INTRODUCTION 


It  is  generally  believed  that  future  improvements  in  the  understanding  of  solid  propellant 
combustion  must  come  from  a  clearer  picture  of  the  chemistry  and  physics  occurring  in  high 
pressure  solid  propellant  flames.  Condensed  phase  reactions  also  need  more  study.  The  structure 
and  reactions  of  solid  propellant  flames  are  not  well  understood  under  rocket  motor  conditions 
(pressures  on  the  order  of  7  MPa  (1000  psia)),  although  enough  is  known  that  serviceable  solid 
propellant  combustion  models  have  been  created  These  models  tend  to  emphasize  the  heat  transfer 
aspects  of  combustion  at  the  expense  of  the  chemistry.  This  allows  much  useful  information  to  be 
derived  from  the  models,  although  many  situations  that  are  directly  chemistry  related,  such  as 
binder  structure  influences  on  burning  behavior,  are  not  well  handled.  To  improve  the 
understanding  of  solid  propellant  combustion,  much  more  detailed  information  is  needed  on 
propellant  flames,  preferably  under  high  pressure  combustion  conditions.  It  is  important  to 
understand  the  combustion  chemistry  since  this  is  the  type  of  information  the  propellant  formulator 
or  chemist  can  use  to  tailor  and/or  improve  the  performance  of  the  propellant. 

The  temperature  profile  in  a  solid  propellant  flame  is  an  important  aspect  of  propellant 
combustion.  Most  state-of-the-art  propellant  combustion  models  (Ref.  1,2)  are  based  upon  a  heat 
transfer  analysis,  with  the  heat  released  in  the  gas  phase  and  condensed  phase  reactions  being 
balanced  against  the  heat  needed  to  bring  the  propellant  from  the  bulk  temperature  to  the  surface 
temperature  plus  the  heat  needed  for  any  endothermic  reactions  such  as  surface  pyrolysis.  In  such 
an  analysis,  the  gas  phase  temperature  profile  would  be  a  key  piece  of  data,  allowing  calculation  of 
the  gas  phase  flame  structure  and  heat  release  (e.g.,  see  Ref.  3).  The  dependence  of  the 
temperature  profile  on  pressure  gives  many  clues  about  the  flame  chemistry,  as  well  as  giving  an 
important  test  of  any  combustion  model.  A  temperature  profile  would  be  a  good  test  of  detailed 
kinetic  models  of  propellant  flames  now  appearing  (Ref.  4). 

Historically,  temperature  profiles  in  propellant  flames  (and  in  the  condensed  phase  as  well) 
have  been  made  with  microthermocouples  (Ref.  3,5,6).  The  thermocouple  beads  (often  <  10  pm) 
are  imbedded  into  the  propellant.  As  the  propellant  bums  down  past  the  thermocouple,  the 
temperature  profile  is  obtained.  Thermocouple  temperature  measurements  have  several  advantages 
and  disadvantages.  Among  the  advantages  are  the  relatively  simple  (inexpensive!)  experimental 
apparatus  needed  and  the  ability  to  measure  the  temperature  in  both  the  condensed  and  gas  phases. 
Among  the  disadvantages  are  the  relatively  limited  temperature  range  (Tmax~  2000  K),  the 
possibility  of  perturbing  the  flame  since  the  thermocouple  is  an  intrusive  probe,  and  the  possible 
problems  of  the  thermocouple  response  time  and  correction  for  heat  loss.  The  environment  of  a 
solid  propellant  flame  is  a  very  difficult  one,  with  temperature  changes  of  -  lO4  K/s  occurring  near 
the  surface.  It  is  this  problem  of  thermocouple  response  time  (resolution)  that  leads  to  an  interest 
in  laser-based  combustion  diagnostics  for  temperature  measurements.  These  techniques  are,  of 
course,  limited  to  the  gas  phase,  although  a  similar  technique  has  been  developed  for  surface 
temperature  measurements  (Ref.  7). 

Laser-based  combustion  diagnostics  such  as  coherent  anti-Stokes  Raman  scattering 
(CARS)  and  laser-induced  fluorescence  (LEF)  are  often  used  to  make  temperature  measurements  in 
flames  (Ref.  8).  CARS  can  be  used  to  measure  the  temperature  of  major  flame  species  such  as  CO 
and  N2,  while  LEF  is  used  to  measure  the  temperature  of  radicals  such  as  CN  and  OH.  Seeding 
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atomic  species  into  flames  to  make  temperature  measurements  using  LEF  is  also  common  (Ref. 
9,10)  (this  type  of  seeding  would  be  relatively  easy  in  a  solid  propellant  flame).  The  LIF  rotational 
temperature  obtained  by  these  techniques  is  that  of  the  ground  electronic  state  of  the  molecule, 
which  is  generally  equal  to  the  translational  temperature  of  the  flame.  Thus  LIF  temperatures  do 
not  yield  the  high  excitation  temperatures  calculated  from  emission  (chemiluminescence)  spectra 
(Ref.  11,12),  which  can  be  well  above  the  adiabatic  flame  temperature.  This  report  describes  an 
attempt  to  apply  OH  LIF  temperature  measurements  to  high  pressure  solid  propellant  flames. 
Successful  OH  LIF  temperature  measurements  have  been  made  in  atmospheric  pressure  propellant 
flames  (Ref.  13).  Successful  CARS  temperature  measurements  on  high  pressure  propellant  flames 
have  recendy  been  reported  (Ref.  14).  Other  spectroscopic  methods  have  been  used  to  measure 
propellant  flame  temperatures,  such  as  sodium  line-reversal  (Ref.  15)  and  infrared  emission  (Ref. 
16).  The  main  limitation  of  these  techniques  is  that  the  temperature  obtained  is  a  Iine-of-sight 
average  of  emission  and  thus  only  moderate  spatial  and  temporal  resolution  is  obtained.  In 
general,  then,  only  thermocouples  have  generated  useable  temperature  measurements  in  solid 
propellant  flames  until  recently. 

LIF  measurements  in  solid  propellant  flames  are  generally  difficult  because  of  the  optical 
thickness  of  the  propellant  flames  (Ref.  12,18).  The  scattering  from  the  panicles  in  the  flames  is 
of  much  higher  intensity  than  LIF  at  the  same  or  nearby  wavelengths,  so  resonance  fluorescence 
(LEF  collected  at  the  excitation  wavelength)  experiments  are  essentially  impossible.  Exciting  and 
collecting  the  fluorescence  in  different  bands  of  a  molecule  is  the  usual  solution  for  this  problem. 
For  example,  one  could  excite  CN  LIF  in  the  0,1  band  and  observe  LIF  in  the  0,0  band  (Ref. 
12,18).  The  OH  LIF  temperature  measurements  described  in  this  paper  were  performed  by 
exciting  LEF  in  the  OH  1,0  band  at  —281  nm  and  observing  the  resulting  LEF  in  the  OH  1,1  band 
from  312-320  nm.  Filters  can  also  be  used  to  aid  in  rejection  of  unwanted  scattered  light  at  the 
laser  wavelength.  The  particles  that  cause  the  scattering  can  also  cause  severe  beam  attenuation 
(Ref.  12,18).  With  certain  propellants,  this  can  occur  at  relatively  low  pressures  (1-5  atm). 
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EQUIPMENT  AND  PROCEDURES 


The  solid  propellant  flames  were  studied  by  burning  propellant  strands  in  a  nitrogen- 
purged,  high  pressure  combustor  (Ref.  17).  This  combustor  is  a  variant  of  the  standard  "window 
bomb"  or  "strand  burner"  in  widespread  use  in  solid  propellant  combustion  studies.  Optical  access 
is  provided  by  eight  3.8-cm  diameter,  2-cm  thick  sapphire  windows  distributed  at  45  degree 
intervals  around  the  combustor.  The  diagnostic  laser  intensity  (and  presumably  any  LIF  signals 
generated)  was  reduced  20%  in  passing  through  one  of  the  windows.  For  laser-diagnostic 
purposes,  a  means  to  keep  the  burning  propellant  surface  at  a  constant  level  is  often  desirable. 
This  is  accomplished  by  a  servo  mechanism  which  elevates  the  surface  of  the  solid  propellant  as  it 
bums  down,  thus  keeping  the  surface  at  a  constant  level  relative  to  the  collection  volume  of  the 
detection  systems.  The  servo  system  employs  a  100  mW  HeNe  laser  passing  across  the  propellant 
surface,  with  a  feedback-control  system  elevating  the  propellant  strand  to  keep  the  fraction  of  the 
beam  blocked  by  the  surface  constant  (Ref.  17).  In  general,  the  system  works  fairly  well  at 
pressures  below  2  MPa  (300  psia),  depending  on  the  propellant,  but  the  increases  in  flame  optical 
thickness  and  propellant  bum  rate  that  accompany  increases  in  pressure  cause  the  system  to  work 
poorly  at  higher  pressures.  An  alternative  is  to  allow  the  propellant  to  bum  down  through  the 
collection  volume,  in  essence  collecting  sequential  data  at  increasing  distances  above  the  propellant 
surface.  All  experiments  were  run  under  nitrogen.  The  propellant  was  burned  in  strands  of  about 
3-4  cm  in  length  and  6  mm  in  diameter,  with  the  sides  of  the  propellant  coated  with  a  thin  layer  of  a 
fluorocarbon  grease  as  an  inhibitor.  Depending  on  the  pressure,  these  strands  burned  from  1  to  20 
seconds.  The  short  bum  time  of  the  strands  did  not  permit  spectral  scanning;  the  detection  system 
was  based  on  the  need  for  rapid  data  collection.  A  premixed  flame  source  constructed  from  6  mm 
tubing  could  be  placed  in  the  combustor  in  place  of  the  strand  for  tuning  of  the  lasers  and  alignment 
of  the  detection  systems.  The  flame  was  usually  run  with  a  N2O/CH4  (molar)  ratio  of 
approximately  2.  Several  propellants  were  examined:  HMX1  (70%  HMX,  17%  TMETN,  13% 
polyester  binder),  HMX2  (80%  HMX,  20%  polyester  binder),  API  (87%  AP,  13%  polybutadiene 
binder),  and  AN1  (72%  AN,  12%  TMETN,  16%  glycidyl  azide  polymer  binder).  The  propellant 
bum  rates  were  HMX1,  1.4  mm/s  at  1.5  MPa  (200  psig);  HMX2,  0.8  mm/s  at  1.5  MPa  (200 
psig);  API,  2.4  mm/s  at  0.5  MPa  (50  psig);  and  AN1, 1.7  mm/s  at  0.8  MPa  (100  psig).  For  more 
information  (see  Ref.  12,18). 

The  OH  LIF  experiments  were  performed  by  pumping  a  Lambda-Physik  FL  2002  dye  laser 
with  308  nm  light  from  a  Lambda-Physik  EMG  201  MSC  excimer  laser.  The  laser  system  had  a 
maximum  repetition  rate  of  80  Hz.  The  dye  laser  was  operated  with  Lambdachrome  LC  5400  dye 
(Coumarin  153),  with  a  nominal  bandwidth  of  0.1  cm'1  after  frequency  doubling  and  a  pulse  width 
of  25  ns.  The  beam  was  focussed  into  the  propellant  flame  with  a  10-cm  focal  length  fused  silica 
lens.  The  resulting  focal  volume  diameter  was  approximately  50  Jim.  The  short  focal  length  lens 
minimized  damage  to  the  sapphire  windows.  After  passing  through  the  combustor  windows,  a 
small  portion  of  the  beam  passed  onto  a  photomultiplier  tube  for  monitoring  beam  transmission 
through  the  combustor  and  propellant  flame.  The  laser  power  was  also  measured  before  the 
combustor  with  a  photo-multiplier  tube  to  monitor  laser  pulse  energy  levels.  This  is  illustrated  in 
Figure  1.  In  experiments  where  the  laser  was  scanned,  the  laser  wavelength  changed  in  0.0004 
nm  steps,  with  the  wavelength  held  at  each  particular  value  for  0.33  sec.  Thus  a  scan  over  0.2  nm 
required  approximately  ten  minutes.  In  these  initial  experiments,  the  data  collection  was  not 
synchronized  with  the  laser  scanning  mechanism,  so  data  collection  occurred  while  the  laser 
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wavelength  was  changing.  This  affected  relatively  few  data  points  so  that  the  lack  of 
synchronization  wasn't  a  problem  in  these  experiments.  In  future  experiments,  the  scanning  and 
data  collection  will  be  synchronized. 

The  LIF  signals  were  collected  at  90  degrees  from  the  laser  beam,  collimated,  and  focussed 
onto  a  spectrometer  entrance  slit  with  1:1  magnification  using  two  10  cm  focal  length  fused  silica 
lenses.  The  image  was  rotated  90  degrees  with  a  mirror  assembly  such  that  the  vertical 
spectrometer  entrance  slits  were  projected  horizontally  onto  the  propellant  flame.  The 
photomultiplier  tube  (PMT)  detection  system  employed  a  RCA  C31034A  tube  attached  to  a  Spex 
1870  1/2  meter  monochromator  (1200  and  300  grooves/mm  gratings,  ~1.5  and  6  nm/mm 
dispersion).  The  collected  LIF  signals  were  focussed  onto  the  slit  of  the  monochromator,  with  the 
vertical  resolution  of  the  system  in  the  flame  thus  being  the  entrance  slit  width  (typically  20-50 
pm).  The  slit  height,  and  thus  the  width  of  the  collection  volume  along  the  path  of  the  laser  beam, 
was  typically  2  mm.  In  earlier  LIF  experiments  (Ref  12,18),  it  was  often  found  necessary  to  add  a 
narrowband  (~10  nm  bandwidth)  filter  to  aid  in  the  separation  of  the  LIF  signals  and  the  scattered 
laser  light  in  the  propellant  flames.  The  PMT  signal  was  processed  by  a  Stanford  Research 
Systems  SR250  gated  integrator,  along  with  the  laser  pulse  energy  and  transmission  signals.  The 
data  collection  was  controlled  by  an  Apple  Macintosh  computer  through  a  Stanford  Research 
Systems  SR245  computer  interface.  Post  processing  was  accomplished  on  a  Macintosh  computer 
using  Kaleidagraph  software.  Propellant  experiments  were  usually  performed  by  setting  the 
spectrometer  on  the  desired  wavelength  and  monitoring  the  LIF  signal  as  a  function  of  time,  as  the 
propellant  burned  down  through  the  collection  volume,  thus  single-pulse  data  acquisition  and  high 
laser  repetition  rates  were  desirable  to  resolve  the  propellant  flame  structure  as  finely  as  possible. 

Some  OH  emission  results  are  also  mentioned  in  this  report  for  comparison  purposes.  In 
these  experiments,  detection  was  performed  with  an  EG&G/PAR  1420  Reticon  diode  array  (700 
elements)  attached  to  a  Spex  1877  Triplemate  spectrometer.  The  Reticon  was  controlled  by  an 
EG&G/PAR  OMA II 1215/1218  detector  controller  and  data  acquisition  system.  The  Triplemate  is 
a  triple  spectrometer  with  an  initial  "filter"  stage  which  allows  high  light  rejection  (Ref.  12,18). 
Three  turret-mounted  gratings  were  available,  with  the  highest  resolution  grating  being  a  2400- 
gr/mm  grating  with  a  dispersion  of  0.0175  nm/pixel.  The  Reticon  was  used  extensively  in  studies 
of  solid  propellant  flame  chemiluminescence,  where  accurate  background  measurements  were 
imperative  due  to  the  presence  of  a  strong  continuous  background,  resulting  from  CO+O 
chemiluminescence  or  emitting  particles.  The  sensitivity  of  the  Reticon  was  usually  too  low  for 
single-pulse  LIF  experiments  (Ref.  12,18). 

The  0.5-m  spectrometer/PMT  system  has  several  difficulties  in  operation  in  propellant 
flames.  Heavily  particulate  laden  flames  cause  Mie  scattering  (or  laser-modulated  particulate 
incandescence  (Ref.  19)  or  rotational  Raman  (Ref.  8))  interference  with  the  LIF  signals,  even  with 
a  filter  in  the  optical  system.  Obviously,  the  lack  of  a  simultaneous  measurement  of  a  non-resonant 
wavelength  as  a  background  correction  is  a  real  limitation  in  this  type  of  system.  Separate 
experiments  were  performed  with  either  the  laser  or  spectrometer  detuned  from  resonance  to  assess 
the  contribution  of  the  non-resonant  scattering  to  the  desired  LIF  signals  (Ref.  12,18). 
Additionally,  the  optical  thickness  of  the  propellant  flames  increases  dramatically  with  pressure 
(higher  soot  density?),  so  higher  pressure  measurements  are  doubly  difficult  (Ref  12,17,18). 
However,  for  several  propellants  with  relatively  non-sooting  flames,  unambiguous  CN  and  OH 
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LIF  signals  were  obtained  up  to  pressures  of  3.5  MPa  using  this  system.  For  very  sooty 
propellants  such  as  AN1,  the  desired  signal  was  difficult  to  obtain  under  most  circumstances.  The 
origin  of  the  soot  is  probably  found  in  the  fuel-rich  stoichiometry  of  the  solid  propellants.  In 
studies  on  a  wide  variety  of  hydrocarbons,  Glassman  found  <{>c's  (critical  equivalence  ratios  above 
which  a  flame  produced  soot)  of  1±0.2  at  flame  temperatures  of  1800-2400  K  (Ref.  46).  The 
stoichiometric  equivalence  ratio  was  based  on  CO  and  H2O.  It  was  also  found  that  (1)  the 
tendency  to  soot  increased  (<j>c  decreased)  for  a  given  stoichiometry  as  the  temperature  decreased, 
(2)  <j)c  decreased  as  the  "number  of  carbon-carbon  bonds"  in  the  fuel  molecule  increased,  and  (3) 
diffusion  flames  and  premixed  flames  show  different  sooting  behavior  for  a  given  fuel.  Using  a 
modification  of  Jain's  method  of  calculating  propellant  stoichiometries  (Ref.  47)  to  account  for  a 
definition  of  stoichiometric  based  on  CO  and  H2O,  the  stoichiometries  of  the  propellants  in  this 
paper  were  calculated  as  4>e  =  1.23,  1.46,  0.95,  and  0.90  for  HMX1,  HMX2,  API,  and  AN1, 
respectively.  <J)  >  1  implies  fuel-rich  stoichiometry,  4>  =1  implies  stoichiometric  (CO  and  H2O), 
and  <j)  <  1  implies  oxidizer-rich  stoichiometry.  All  the  propellants  are  within  the  range  of 
stoichiometries  at  which  flames  have  been  found  to  soot,  and  the  lower  temperature  of  the  AN  1 
flame  would  help  explain  the  increased  sooting  tendency  of  this  flame  (as  evidenced  by  the  optical 
thickness  of  the  flame)  compared  to  the  AP  flame.  If  the  sooting  tendency  of  propellant  flames 
could  be  related  to  the  fuel  decomposition  mechanisms  (as  has  been  done  for  hydrocarbons  (Ref. 
46)),  some  interesting  results  would  undoubtedly  emerge.  For  laser-diagnostic  purposes, 
however,  the  soot  is  a  large  hindrance  to  obtaining  species  and  temperature  data. 
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Figure  1.  Schematic  (top-view)  of  apparatus. 
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TECHNIQUE  APPLICATION 


OH  is  found  in  all  flames  containing  hydrogen  and  oxygen  and  is  often  found  in  high 
concentrations.  For  example,  the  equilibrium  mole  fraction  of  OH  in  a  stoichiometric  CH4/N2O 
(1:4)  flame  at  1  atm  is  0.04.  For  comparison,  the  equilibrium  OH  mole  fractions  and  adiabatic 
flame  temperatures  (K)  of  the  four  propellants  at  3.5  MPa  (500  psig)  are  0.0142/2970  (API), 
0.00066/2620  (HMX1),  8X10-8/2050  (HMX2),  and  0.0001/2140  (AN1).  The  equilibrium  OH 
concentration  is  not  the  absolute  concentration  expected,  but  is  an  indication  of  the  relative 
stoichiometries  in  the  various  flames.  In  addition  to  occurring  at  high  concentration,  OH  is  a  good 
temperature  probe  since  it  is  usually  distributed  throughout  a  flame  (Ref.  20,21).  This  was  also 
found  to  be  the  case  in  several  propellant  flames  (Ref.  12).  In  contrast,  radicals  such  as  CN  and 
NH  are  confined  to  certain  propellant  reaction  zones,  typically  very  near  the  surface,  and  thus  are 
not  as  versatile  as  OH  as  a  temperature  probe.  Finally,  OH  is  probably  the  best  understood  radical 
in  a  flame  environment,  with  many  studies  made  of  the  spectroscopy  and  laser-excitation  dynamics 
of  OH  (e.g.,  Ref.  22-37).  For  certain  propellants  (notably  double-base  propellants  and  low- 
energv  nitramine  propellants  which  have  "dark"  zones  where  one  would  expect  little  OH)  and 
certain  temperature  regimes  (low  temperatures  especially),  OH  is  probably  not  the  best  molecule 
for .  .IF  studies.  A  good  possibility  for  these  propellants  and  conditions  is  NO.  The  possibility  of 
using  NO  LIF  as  a  temperature  probe  in  propellant  flames  is  under  investigation  at  NWC.  NO  has 
been  used  as  a  flame  thermometer  (e.g.,  Ref.  38).  NO  has  several  unique  problems,  notably  the 
high  concentrations  found  in  propellant  flames,  causing  beam  attenuation/absorption  problems,  and 
the  relative  (compared  to  OH)  lack  of  spectroscopic  information. 

OH  LIF  measurements  are  performed  by  exciting  transitions  in  the  OH  A2Z-X2n  system 
(Ref.  39).  For  applications  where  the  separation  of  the  desired  fluorescence  from  scattered  laser 
light  is  difficult,  such  as  2-D  imaging  (Ref.  13,40)  or  propellant  flame  probing  (Ref.  13,18),  the 
LIF  excitation  and  detection  are  performed  in  different  bands  separated  by  20-30  nm.  For 
example,  the  measurement  might  be  excite  0,1 -detect  0,0  or  excite  1,0-detect  1,1.  The  spectral 
separation  of  the  bands  enhances  the  separation  of  the  desired  LEF  from  the  undesired  scattering  at 
the  laser  wavelength.  For  these  experiments,  the  second  measurement  method  (excite  1,0-detect 
1,1)  is  used.  For  example,  the  LEF  spectra  obtained  in  the  outer  cone  of  the  1  atm  1:2  CH4/N2O 
calibration  flame  from  the  excitation  of  the  R^  5  line  of  the  OH  1,0  band  are  shown  in  Figures  2 
(0,0  and  1,1  band)  and  3  (1,0  band).  A  schematic  of  the  LIF  is  shown  in  Figure  4.  In  the  case 
where  no  rotational  redistribution  occurs  in  the  laser-populated  upper  state,  the  only  lines  present  in 
the  spectra  would  be  the  Rt  5,  Qj  6,  and  7  lines  in  the  1,0  and  1,1  bands  (ignoring  weak 
satellite  transitions).  As  can  be  seen  in  Figures  2  and  3,  considerable  rotational  redistribution  has 
occurred  in  the  upper  state  prior  to  the  molecule  radiating.  Vibrational  redistribution  has  also 
occurred,  as  shown  by  the  presence  of  0,0  band  LIF.  The  excite  1,0-detect  1,1  scheme  has  the 
added  advantage  of  minimizing  the  self-absorption  problems  that  occur  when  the  0,0  band  is 
involved  (Ref.  34,36). 

This  rotational  redistribution  makes  calculations  of  OH  concentration  more  complicated.  In 
a  well-characterized  flame  environment,  detailed  calculations  involving  the  dynamics  of  the  LEF 
have  been  made  (Ref.  8).  Saturating  the  OH  LIF  simplifies  these  calculations  (Ref.  8).  For 
temperature  measurements,  detailed  calculations  are  typically  avoided  by  using  broadband 
collection;  for  example,  collecting  all  the  fluorescence  from  the  1,1  band  in  Figure  2.  A  bandpass 
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of  ~20  nm  is  required  (Ref.  28).  This  was  obtained  in  our  experiments  by  replacing  the  1200- 
gr/mm  grating  used  in  earlier  LIF  work  (Ref.  12,18)  with  a  300-gr/mm  grating,  which  gave  a 
bandpass  of  ~20  nm  when  the  rear  slit  of  the  spectrometer  was  opened  to  3  mm.  This  is  shown  in 
Figure  5,  where  the  bandpass  is  measured  by  scanning  the  spectrometer  over  the  laser  line.  The 
bandpass  was  centered  at  315  nm  for  the  OH  LIF  temperature  measurements,  as  illustrated  in 
Figure  6.  Temperatures  are  obtained  in  this  type  of  experiment  by  scanning  the  laser  across  many 
rotational  lines  in  the  linear  or  saturated  fluorescence  regime  and  collecting  the  resulting 
fluorescence  (an  "excitation"  scan  (Ref.  8)).  The  collected  fluorescence  is  proportional  to  the 
ground  state  population  of  the  given  transition,  the  absorption  coefficient  (transition  probability), 
and  the  degeneracy  of  the  ground  state  rotational  level.  If  one  assumes  a  Boltzmann  distribution 
(rotational  population  proportional  to  exp  (-rotational  energy/kT)),  then  the  temperature  can  be 
found.  Thus,  Ifi'-g"B(v,J)exp(-hcEr/kT),  where  If]  is  the  fluorescent  intensity,  g"  is  the 
degeneracy  of  the  ground  state  rotational  level  (g"=2J"+l),  B(vJ)  is  the  Einstein  coefficient  for 
absorption  (Ref.  41),  and  Er  is  the  rotational  energy  of  the  ground  state  rotational  level  {found 
from  tables  (Ref.  39,48)  or  approximated  as  ~BVJ”(J"+1)  (ref.  8)).  The  data  used  are  listed  in 
Appendix  1.  Thus  a  plot  of  In  If]/g"B(v  J)  vs  Ef  will  have  a  slope  of  -hc/kT.  These  equations  are 
valid  for  the  linear  fluorescence  regime,  where  LIF  signal  is  proportional  to  laser  power .  Under 
saturation  conditions,  the  equations  change  somewhat  (Ref.  35),  with  -hc/kT  being  the  slope  of  a 
plot  of  In  [Ifl(l/g"+l/g')l  vs  Er. 

An  example  of  an  excitation  scan  is  shown  in  Figure  7.  Assignment  of  lines  was  simple 
and  accurate,  as  shown  in  Figure  8.  The  estimated  spectral  irradiance  for  the  experiment  of  Figure 
7  is  approximately  8X106  W/cm2-cm-1,  much  larger  than  the  saturation  estimate  of  1X106  W/cm2- 
cnr1  for  atmospheric  pressure  flames  (Ref.  8).  The  irradiance  was  calculated  using  a  bandwidth  of 
0.1  cm*1  (dye  laser  before  frequency  doubling  is  0.2  cm'1),  a  pulse  duration  of  25  ns,  and  a 
focussed  beam  diameter  of  50  Jim.  The  focussed  beam  diameter  was  estimated  as  d=f0,  with  f 
(lens  focal  length)=10  cm  and  ©  (beam  divergence)=0.5  mrad.  When  the  data  was  reduced  as 
described  assuming  complete  saturation,  it  was  found  that  the  various  branches  did  not  fall  on  a 
single  line.  This  is  shown  for  the  excitation  scan  of  Figure  7  in  Figure  9,  yielding  best  fit 
temperatures  of  2930,  1584,  and  2638  K  for  the  R],  Rj',  and  R2  branches.  If  the  same  data  is 
reduced  by  assuming  non-saturation,  a  similar  type  of  Boltzmann  plot  is  obtained;  as  shown  in 
Figure  10,  best  fit  temperatures  of  2373,  3560,  and  2288  K  are  obtained  for  the  Rj,  Rj',  and  R2 
branches.  In  both  cases,  the  main  branches  (R]  and  Rj)  fall  on  essentially  the  same  curve  (as 
shown  in  Figure  10),  while  the  satellite  R]'  branch  gives  a  substantially  different  temperature.  A 
consistent  exception  to  this  linear  behavior  was  seen  for  those  lines  with  low  quantum  numbers, 
such  as  R]  1  and  Rj  2,  which  were  ignored  in  the  temperature  fits  discussed  above.  Non-linearity 
of  Boltzmann  plots  at  low  rotational  energy  is  relatively  common  (e.g..  Ref.  49). 

When  a  similar  experiment  was  performed  with  a  higher  laser  energy  (14  pJ/pulse),  the 
results  shown  in  Figures  12  and  13  were  obtained.  The  OH  LIF  temperatures  calculated  were 
2763,  1535,  and  2534  K  (saturated  assumption)  and  1905,  3695,  and  1954  K  (non-saturated 
assumption)  for  the  Rj,  Rj',  and  R2  branches.  It  is  apparent  that  the  OH  LIF  is  still  partially 
saturated,  even  with  the  substantially  higher  laser  energy.  In  fact,  saturation  was  not  evident  in  an 
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Figure  4.  OH  LIF  schematic. 


LIF-vs-laser  power  plot  until  laser  powers  on  the  order  of  100  pJ/pulse  was  reached,  as  shown  in 
Figure  14.  One  possible  reason  that  the  experimental  data  is  partially  saturated,  even  though  the 
irradiance  is  greater  than  the  saturation  estimate,  is  that  saturation  is  sensitive  to  the  detection  gate 
used  (both  width  and  delay)  and  other  parameters  such  as  inclusion  of  the  "wings"  of  the  laser 
beam  (Ref.  8,35).  Using  a  method  for  reducing  partially  saturated  data  (Ref.  45),  the  temperature 
for  the  excitation  scan  of  Figure  7  is  260Q±80  K,  with  the  Rj,  Rj’  and  R2  branches  falling  on  the 
same  curve  (calculation  performed  by  Steve  Zabamick,  UDRI).  The  adiabatic  flame  temperature  of 
the  calibration  flame  is  250G±400  K  (the  uncertainty  is  due  to  uncertainty  in  the  gas  flow  rates  and 
the  flame  stoichiometry).  Thus,  the  corrected  temperature  is  apparently  reasonable.  This  data  may 
be  compared  with  that  of  Anderson,  et  al  (Ref.  22),  where  all  branches  appeared  to  fit  on  one  line 
in  an  In  Ifl/g"B(v,J)  vs  E,.  plot  (in  an  excite  1,1-detect  1,1  experiment),  thus  indicating  non¬ 
saturation.  The  spectral  irradiance  was  approximately  3X106  W/cm^cnr1  (1  (is  pulse,  0.34  cm*1 
bandwidth,  200  pm  beam  focussed  diameter,  0.3  mJ/pulse,  1 1.4-cm  focal  length  focussing  lens). 
This  discussion  indicates  the  importance  of  saturation  corrections  in  obtaining  accurate  OH  LIF 
temperature  measurements,  especially  when  the  satellite  bands  are  used. 

Although  the  calibration  flame  used  in  these  experiments  was  not  the  well-characterized 
flame  needed  for  detailed  spectroscopy  studies,  it  appeared  valuable  to  undertake  some  background 
experiments  to  examine  the  effect  of  various  parameters,  such  as  gate  width  and  timing,  on  LIF 
temperature  measurements.  It  is  anticipated  that  detailed  studies  will  be  performed  in  the  high 
pressure  flame  facility  under  construction  in  this  laboratory.  In  addition,  the  effect  of  a 
narrowband  filter  (Oriel  56510,  313  nm  center  bandpass,  11-nm  FWHM  bandwidth,  25%  peak 
transmission)  on  the  temperature  measurements  was  examined  for  use  during  the  propellant  flame 
OH  LIF  temperature  measurements. 
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Wavelength,  nm 


ndpass  of  "broadband"  detection  system.  Measured  by  scanning  spectrometer 
laser  line  (281.13  nm),  300  gr/mm  grating,  50-^irn/3-mm  slits,  no  flame. 
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CH4/N2O  flame,  1  atm  (air),  outer  cone,  0.4  fiJ/pulse  in  flame,  10  Hz,  3 
pulse  average,  30  pm/3  mm  slits,  300  gr/mm  grating,  43000  data  points. 
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Figure  8.  Comparison  of  measured  line  positions  of  Figure  ?  with  literature  yah 

(Ref.  39). 
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Figure  9.  Boltzmann-tvpe  of  plot  of  data  from  Figure  7.  assuming  saturation. 


Figure  10. 


;igure  11. 


* 

-j _ 1 _ 1 _ 

l  .  1 _ i _ 1  »  -  1 _ i _ 1 _ i _ 1 _ i _ 

1 1/88-43to49 

0°D  □ 

-5.8  - 

Q> 

□  D 

9>  O 

□  □ 

-6.2  - 

o 

-6.6  - 

_i _ 

□ 

□ 

6 

-7  - 

s 

-7.4  - 

o  R1 

°  $ 

□  R 1 

9> 

-7.8  - 

o  R2 

o  ~ 

-8.2  - 

— 

- - 1 

-i - -  r  — i  —  ■  t  ■  i  ~  »  ~ 

0.000  500.0  1000  1500  2000  2500  3000  3500  4000 

Rotational  energy,  i/cm 

Boltzmann-tvpe  plot  of  data  from  Figure  7.  assuming  non-saturation. 


Rotational  energy,  1/cm 


Temperature  fits  to  data  from  Figure  10.  Data  with  N"= 1,2  eliminated. 
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Figure  14.  Saturation  assessment.  Broadband  315  nm  detection,  281.915  nm  excitation 
(OH  1,0  Ch  1  +  R2  3).  50  |im/3  mm  slits,  CH4/N2O  flame,  outer  cone,  1  atm 

(air). 


The  effects  of  varying  the  excitation  line  on  the  structure  of  the  LIF  spectrum  are  shown  in 
Figures  15-17.  It  can  be  seen  that  all  the  strongly  excited  lines  arc  typically  within  the  bandpass  of 
Figure  6.  To  be  outside  the  bandpass  (>324  nm),  the  observed  1,1  band  lines  would  have  to  be 
Pj,  P2  (N">17)  or  Qj,  Q2  (N">24).  The  structure  of  this  LIF  spectrum  is,  of  course,  radically 
different  from  the  OH  emission  spectrum  from  the  flame,  an  example  of  which  is  shown  in  Figure 
18.  For  comparison  purposes,  the  LIF  spectra  in  the  0,0,  1,0,  and  0,1  bands  resulting  from  OH 
0,0  excitation  are  shown  in  Figures  19-21.  The  narrowband  filter  has  an  important  effect  on  the 
structure  of  the  LEF  reaching  the  detector.  This  is  shown  in  Figures  22  and  23,  where  the  emission 
spectrum  from  the  flame  is  shown  with  and  without  the  filter  present.  Note  that  the  filter  attenuates 
the  R  branches  of  the  0,0  band  severely.  The  filter  also  affects  the  1,1  band  slightly,  as  can  be 
seen  by  comparing  Figures  15  and  24.  The  effect  the  filter  has  on  LIF  temperature  measurements 
is  not  immediately  obvious,  so  the  filter  was  used  during  some  of  the  excitation  scans  described  in 
the  next  paragraph  to  determine  the  influence  of  the  filter. 
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Figure  15.  OH  LIF  spectrum.  1.0  Rt  5  excitation.  Conditions  similar  to  Figure  2  with 

300-gr/mm  grating. 
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Figure  16.  OH  LIF  spectrum.  1.0  Rt  2  excitation.  Conditions  same  as  Figure  15. 
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Figure  17.  OH  LIF  spectrum.  l.QR^  8  excitation.  Conditions  same  as  Figure  15. 


Figure  18.  OH  emission  spectrum.  CH4/N2O  flame,  Reticon  detector,  700  pixels,  0.0175 

nm/pixel  dispersion,  500/200  |xm  slits,  0.2  sec  exp.  time. 
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Figure  19.  OH  Q.O  UF  spectrum.  0.0  R,  8+R}  10  excitation  (306.36  nm).  CH4/N2O 
flame,  1  atm  (air),  outer  cone,  0.03  mJ/pulse  in  flame,  50  pm  slits,  1200- 

gr/mm  grating. 


Figure  20. 


l.  Same  conditions  as  Figure  19. 
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Figure  21.  OH  1.0  LIF  spectrum.  Same  conditions  as  Figure  19. 


Figure  22.  OH  emission  spectrum.  Same  conditions  as  Figure  18  except  1200-gr/mm 

grating  (0.035  nm/pixel  dispersion). 
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OH  emission  intensity 


To  examine  the  effect  of  various  detection  parameters  on  OH  LIF  temperature 
measurements,  excitation  scans  of  the  length  of  that  in  Figure  7  were  too  time-consuming.  For  this 
reason  a  shorter  scan  over  the  OH  1,0  Rj  and  Rj'  branches  from  approximately  281.1  to  281.35 
nm  was  employed.  An  example  of  such  an  excitation  scan  and  the  resulting  Boltzmann  plot  for  the 
Rj  branch  (Rj  2  -Rj  8)  is  shown  in  Figures  25  and  26.  Unfortunately,  this  truncated  scan  gave 
temperatures  that  were  consistently  lower  than  those  of  the  more  complete  scan  of  Figure  7.  For 
example,  using  only  the  Rj  2  to  Rj  8  lines  of  Figures  7  and  10  gave  a  best-fit  temperature  of  1579 
K,  while  using  Rj  3  to  Rj  8  gave  a  best-fit  temperature  of  1749  K  (unsaturated  assumption).  It 
appears  best  to  avoid  low  J"  lines  in  this  type  of  measurement  altogether.  A  more  complete 
analysis  than  the  one  below  would  use  a  larger  range  of  lines  to  obtain  more  accurate  data.  For  the 
purposes  of  this  study,  a  limited  scan  was  deemed  to  be  appropriate.  For  better  accuracy,  the  data 
should  be  corrected  for  partial  saturation  also.  For  the  limited  purposes  of  this  study,  the  data  were 
not  corrected.  Instead,  the  OH  LIF  excitation  scans  were  reduced  for  both  the  saturated  and 
unsaturated  limiting  cases,  and  the  effects  of  detection  gate  width  and  timing  were  assessed  by  the 
changes  in  the  limiting  case  behavior. 

The  temporal  location  of  the  detection  gate  (relative  to  the  laser  pulse)  should  have  an 
important  influence  on  the  characteristics  of  the  OH  LIF  excitation  scans  and  thus  could  have  an 
important  influence  on  the  OH  LIF  temperature  measurements.  In  the  results  discussed  thus  far,  a 
30-ns  detection  gate  was  used,  with  the  timing  adjusted  so  that  the  LEF  intensity  was  a  maximum. 
This  occurred  at  a  gate  delay  (relative  to  a  trigger  pulse  from  the  laser  power  supply)  of  1.20  ps. 
For  a  10-ns  gate  width,  the  maximum  OH  LIF  intensity  occurred  at  essentially  the  same  delay 
setting.  For  a  10-ns  gate  width,  changing  this  delay  setting  by  ±0.01  ps  (10  ns)  reduced  the  LIF 
intensity  by  50%.  The  effect  of  this  change  on  OH  LIF  excitation  scans  is  shown  in  Figures  27 
and  28.  With  the  earlier  detection  gate  (Figure  27),  the  excitation  scan  approached  the  unsaturated 
limiting  case  more  closely  than  the  baseline  excitation  scan  (Figure  26),  as  shown  by  the  closeness 
of  the  unsaturated  Rj  and  Rj'  branches  on  the  Boltzmann  plot.  With  the  later  detection  gate 
(Figure  28),  the  Boltzmann  plots  show  more  saturation.  The  trend  is  steady;  as  shown  in  Figure 
29,  a  summary  of  Figures  26-28.  With  the  short  scans  employed  (N”=3-8),  eliminating  one  of  the 
points  on  the  Boltzmann  plots  often  changed  the  best  fit  temperature  by  hundreds  of  degrees. 
Thus,  the  accuracy  of  these  scans  is  too  low  to  draw  definitive  conclusions  from  this  data.  Within 
the  accuracy  of  this  method,  the  excitation  scans  for  the  10  and  30-ns  baseline  gate  widths  gave 
identical  temperatures.  Changing  the  center  of  the  detection  bandpass  to  312.5  and  317.5  nm 
(from  315  nm)  also  had  a  negligible  effect  on  the  rotational  temperature  found  from  the  excitation 
scans.  In  early  experiments  with  a  3-nm  detection  bandpass,  a  strong  effect  of  the  location  of  the 
detection  bandpass  (314-316  nm)  was  found.  With  the  20-nm  bandpass  used  for  most  of  the  OH 
LIF  temperature  measurements,  this  effect  was  not  seen.  In  most  cases,  the  addition  of  the  313- 
nm  filter  also  had  a  negligible  effect  on  the  excitation  scan  temperature.  A  notable  exception  was 
the  10-ns  "early"  gate,  where  the  addition  of  the  filter  had  a  large  effect  on  the  temperatures,  as  can 
be  seen  by  comparing  Figures  27  and  30.  The  baseline  10  and  30-ns  cases  showed  little  effect  of 
adding  the  filter.  Thus,  under  certain  conditions,  a  narrowband  filter  (or  a  narrower  detection 
bandpass),  may  be  used  if  care  is  taken  to  assure  that  the  narrower  bandpass  is  not  affecting  the 
temperature  measurement.  The  main  conclusion  to  be  drawn  from  this  work  on  the  calibration 
flame  is  that  careful  attention  must  be  paid  to  the  saturation  characteristics  of  the  system  to  obtain 
reliable  temperature  measurements. 
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This  type  of  excitation  scan  requires  several  minutes  as  a  minimum  to  obtain  accurate  data. 
In  transient  experiments  where  the  necessary  time  for  an  excitation  scan  is  not  available,  typically 
two  experiments  are  performed  with  two  excitation  lines  (denoted  by  subscript  1  and  2);  such  that 
kT/hc=(Erl-Er2)/[(Ifi/g"B)i/(Ifl/g"B)2]  in  the  unsaturated  case  and  kT/hc=(Erp 
Er2)/{  i/CIflC l/g”-»- l/g')]2)  f°r  ^e  saturated  case.  In  a  good  example  of  the  two-line 

technique,  Parr  &  Parr  at  NWC  obtained  two-dimensional  images  of  OH  LIF  Tr  in  a  laser-driven 
atmospheric  pressure  propellant  flame  by  exciting  two  lines  in  separate  experiments  and  dividing 
the  resulting  digital  images  (Ref.  13).  Of  course,  any  differences  between  the  two  experiments 
would  lead  to  errors  in  the  temperature  profile. 

The  two-line  procedure  involves  several  assumptions,  aside  from  questions  of  saturation. 
Although  broadband  collection  assures  collection  of  the  LIF  from  an  entire  band  when  done 
properly  (Ref.  28),  if  the  various  upper  rotational  levels  have  different  quenching  behavior 
dependent  on  rotational  level,  the  measurements  would  be  biased  and  give  incorrect  temperatures. 
There  is  some  evidence  that  quenching  does  vary  with  rotational  level  in  some  molecules,  including 
OH  (Ref.  32,42).  To  correct  for  this  problem,  Cattolica  (Ref.  34)  uses  the  two-line  excitation 
method,  but  excites  the  same  upper  state  from  two  different  lower  levels  (e.g.,  Q21  9  and  S21  7); 
thus  the  quenching  behavior  would  be  the  same  for  the  two  lines.  This  approach  avoids  the 
quenching  problem.  However,  it  can  have  temperature  sensitivity  problems  since  the  two  levels 
are  of  relatively  close  rotational  energy;  the  maximum  sensitivity  is  obtained  by  means  of  relatively 
weak  satellite  transitions  (Ref.  34).  This  same  scheme  is  employed  when  making  temperature 
measurements  using  seeded  species  such  as  T1  (Ref.  9,10).  Most  of  the  verification  of  the 
excitation  scan  technique  has  been  done  at  low  pressures,  so  there  is  some  question  as  to  the 
application  of  this  technique  to  high  pressure  flame  diagnostics.  The  effects  of  pressure  on  this 
technique  merits  further  study. 
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Figure  25.  OH  LEF  excitation  scan  with  broadband  collection  centered  at  3 15  nm.  Same 

conditions  as  Figure  7  except  0.002  mJ/pulse  laser  energy  in  flame  and  10  ns  gate 

width. 
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Figure  26. 


Figure  27  except  filter  added  to  collection  optics. 


PROPELLANT  FLAME  EXPERIMENTS 


For  the  propellant  flame  measurements  described  in  this  paper,  it  was  decided  that  the 
greater  temperature  sensitivity  obtainable  by  exciting  strong  transitions  with  differing  upper  levels 
compensated  for  the  errors  involved  in  neglecting  quenching  differences.  The  two  lines  selected 
for  the  propellant  experiments  were  Ri  4  and  Rj  10.  The  temperature  response  of  these  two  lines 
is  illustrated  in  Figure  31.  The  propellant  experiments  were  performed  by  allowing  the  propellant 
to  bum  down  through  the  laser  beam,  as  idealized  in  Figure  32.  This  is  the  same  as  keeping  the 
propellant  surface  at  a  constant  height  and  moving  the  laser  beam  and  collection  optics  up  through 
the  flame  (but  much  simpler).  To  perform  the  temperature  measurement,  average  OH  LIF  profiles 
from  R1  4  and  R1  10  excitation  were  obtained  from  multiple  experiments;  the  resulting  profiles 
were  then  ratioed  and  converted  to  a  temperature  profile  point-by-point  using  Fig.- 31.  Several 
things  should  be  noted  in  Figure  32.  For  a  given  propellant  bum  rate,  the  number  of  data  points  in 
a  given  flame  region  is  determined  by  the  laser  repetition  rate.  For  example,  for  a  propellant  bum 
rate  of  2  mm/s,  a  40  Hz  laser  repetition  rate  implies  that  the  propellant  surface  has  regressed  50  pm 
between  data  points.  Thus  a  premium  is  placed  on  high  laser  repetition  rates  in  such  a  "bum- 
through"  experiment.  This  type  of  experiment  is  procedurally  identical  to  a  propellant  micro¬ 
thermocouple  measurement,  except  that  no  data  is  collected  in  the  condensed  phase.  Also,  there  is 
a  region  near  the  propellant  surface  where  there  is  only  partial  transmission  of  the  laser  beam. 
From  a  consideration  of  the  geometry  of  the  experiment,  this  region  should  be  100-200  pm  in 
thickness  for  a  perfectly  flat  propellant  surface.  A  surface  that  is  level  but  "rough",  such  as  the 
surface  of  an  AP  composite  propellant,  would  increase  this  interaction  zone  thickness.  A  non-level 
surface  would  compromise  the  data  further.  As  a  result,  LIF  measurements  inside  this  zone  are 
rather  suspect.  LIF  measurements  in  this  region  could  perhaps  be  normalized  by  the  transmission, 
although  in  optically  thick  flames  this  procedure  does  not  work  well. 


Figure  31.  Temperature  sensitivity  of  OH  1.0  Rt  4  and  RT,  lOexcitation. 
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Figure  32.  Schematic  of  OH  LIF  "bumdown"  experiments. 


Once  the  LIF  profiles  from  Ri  4  and  Rj  10  excitation  are  collected,  a  distance  scale  (above 
the  propellant  surface)  is  assigned;  this  procedure  is  prone  to  some  error  because  the  precise 
location  of  the  surface  is  occasionally  not  obvious.  The  profiles  can  then  be  normalized  as 
described  in  the  previous  section  and  divided,  producing  a  temperature  profile.  For  example,  the 
"OH  LIF"  profile  in  the  flame  above  an  87%  API  propellant  (Ri  4  excitation)  is  shown  in  Figure 
33  for  two  experiments.  Such  profiles  were  obtained  also  in  HMX1  propellant  flame  at  1.5  MPa 
(200  psig)  and  in  the  AN1  propellant  flame  at  0.8  MPa  (100  psig),  as  illustrated  in  Figures  35  and 
36.  The  laser  attenuation  in  the  AN1  flames  was  so  large  that  temperature  measurements  in  this 
propellant  flame  were  abandoned.  The  attenuation  can  be  seen  by  comparing  the  transmission 
scales  of  Figures  1 1  and  12,  taken  with  the  same  laser  power.  Even  in  the  relatively  "clean”  AP 
and  HMX  flames,  the  laser  attenuation  was  often  >50%  at  pressures  above  1-2  MPa  (150-300 
psig)  (Ref.  12,17,18).  The  attenuation  in  the  HMX2  propellant  flame  was  also  very  large.  5-10 
profiles  were  collected  from  both  Rj  4  and  Rj  10  excitation  in  the  API  and  HMX1  propellant 
flames  over  a  period  of  several  days.  The  laser  wavelength  was  periodically  checked  with  the 
calibration  flame  to  ensure  that  the  laser  was  tuned  to  the  correct  line.  The  results  were  relatively 
inconsistent,  with  large  variations  (factor  of  >2)  in  the  peak  "OH  LIF"  intensity,  although  most  of 
the  profiles  were  superficially  similar.  The  variations  from  experiment  to  experiment  can  be  seen 
in  Figure  33.  When  the  laser  was  detuned  from  resonance  (no  OH  LIF  seen  in  the  calibration 
CH4/N2O  flame),  and  an  identical  "bum-through"  experiments  were  performed,  the  results  shown 
in  Figures  37  and  38  were  obtained.  Obviously,  there  is  a  problem.  The  wide  detection  bandpass 
which  ,he  temperature  measurement  requires  was  allowing  large  amounts  of  spurious  laser 
scattering  from  the  flame  to  contaminate  (actually,  overwhelm  is  probably  more  accurate)  the  OH 
LIF  signal,  leading  to  poor  results.  This  result  is  not  surprising  when  one  considers  that  spurious 
scattering  was  a  problem  in  CN  LIF  experiments  where  the  collection  bandwidth  (421.5  nm 
excitation,  388.3  nm  detection)  was  ~0.05  nm  (Ref.  12,18).  Maintaining  a  wide  detection 


27 


bandpass  while  maintaining  a  high  level  of  light  rejection  is  the  major  stumbling  block  to  the 
application  of  this  technique  to  high  pressure  solid  propellant  flames.  One  possibility  would  be  to 
use  a  double  or  triple  spectrometer  to  increase  the  light  rejection.  Another  possibility  might  be  to 
use  a  narrower  bandpass,  carefully  calibrated  (in  the  previously  mentioned  high  pressure  flame)  to 
allow  correction  of  the  narrowband  results  to  yield  accurate  temperatures. 

Additional  propellant  experiments  were  performed  with  the  previously  described 
interference  filter  (313  nm,  1 1-nm  bandpass),  in  the  collection  path  to  aid  in  light  rejection.  As  can 
be  seen  by  comparing  Figure  38  (Rj  10  excitation)  and  39  (Rj  5  excitation)  with  Figure  40  (laser 
detuned  from  resonance),  the  addition  of  the  filter  has  enhanced  the  light  rejection  enough  that  OH 
LIF  is  apparent  above  1  him  in  the  AP  propellant.  Between  the  surface  and  1  mm,  spurious  (non- 
OH  LIF)  scattering  is  still  dominant.  Although  the  results  are  unreliable,  rough  temperature 
estimates  can  be  made  from  2-4  mm  above  the  surface  from  Figures  38  and  39,  yielding  1700- 
1800  K.  Considering  the  biases  discussed  in  the  previous  section,  the  "corrected"  (actual) 
temperature  could  be  ~500  K  higher.  At  a  distance  from  the  surface  of  2-4  mm,  it  seems 
reasonable  to  expect  that  the  flame  temperature  should  be  near  the  adiabatic  flame  temperature  of 
the  propellant  (~2900  K).  Although  the  data  is  suspect,  it  suggests  that  the  flame  temperature  in 
the  API  propellant  flame  is  well  below  the  adiabatic  flame  temperature.  This  discussion  has 
assumed  that  the  API  propellant  flame  is  homogeneous  on  the  scale  of  the  collection  volume, 
which  can  be  approximated  as  a  slab  of  100  |im  in  thickness  and  depth  and  2  mm  in  length.  On 
this  scale,  the  heterogeneous  nature  of  the  propellant  may  persist  into  the  flame.  This  data  can  be 
compared  to  the  temperatures  inferred  from  OH  emission  spectra  in  this  same  propellant  flame  by 
Campbell  (Ref.  1 1).  He  found  that  the  rotational  temperature  of  the  excited  electronic  (A2£)  state 
of  OH  in  the  API  propellant  flame  at  0.8  MPa  at  ~2  mm  above  the  surface  was  290Q±~200  K, 
very  close  to  the  adiabatic  flame  temperature.  As  pressure  increased,  the  emission  temperatures 
also  increased.  His  calculations  were  based  on  comparisons  of  computer-generated  emission 
spectra  and  actual  emission  spectra  in  the  propellant  flame.  The  experimental  spectra  were 
collected  with  the  Reticon  diode  array  detector  as  in  Fig.  18.  The  emission  temperatures  were 
much  more  reproducible  than  the  LIF  temperatures,  although  the  interpretation  of  the  emission 
temperatures  was  made  more  difficult  by  possible  self-absorption  effects  and  non-equilibrium 
rotational  distributions  (Ref.  11).  Also,  the  temperature  of  the  OH  A2Z  state  may  not  be  the  same 
as  that  of  the  much  more  populous  ground  (X^IT)  state.  LIF  experiments  measure  the  temperature 
of  the  X2II  state  (Ref.  8).  It  remains  to  be  proven  whether  or  not  the  OH  A2X  temperature  is  the 
same  as  the  translational  temperature  in  the  flame,  as  is  often  assumed  (Ref.  53,54). 

A  similar  series  of  experiments  was  performed  in  the  HMX1  propellant  flame  at  1.5  MPa 
(200  psig),  as  shown  in  Figures  41-43.  In  this  propellant  flame,  the  OH  LIF  signal  is 
indistinguishable  from  the  scattering.  The  reproducibility  of  the  data  is  also  poor,  as  can  be  seen  in 
Figures  41-43.  The  lack  of  OH  LIF  signal  is  not  surprising,  since  the  equilibrium  OH 
concentration  in  the  HMX1  flame  at  the  adiabatic  flame  temperature  is  about  20  times  smaller  than 
that  in  the  API  propellant  flame.  Thus,  the  HMX1  flame  is  a  more  difficult  subject  for  OH  LIF 
measurements.  Successful  temperature  measurements  have  been  recently  reported  in  a  HMX1 
propellant  flame  using  CARS  (Ref.  14).  The  temperature  several  millimeters  from  the  surface  was 
found  to  be  near  the  adiabatic  flame  temperature  (2600  K). 
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Figure  33.  "OH  LIF"  profile  in  propellant  flame.  API  propellant,  0.5  MPa  (50  psig),  0.015 
mJ/pulse  in  flame,  40  Hz,  Ri  10  excitation  (281.537  nm),  50  pm/3  mm  slits. 
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Figure  34.  "OH  LIF”  profile  in  propellant  flame.  HMX1  propellant,  1.5  MPa  (200  psig),  0.15 
mJ/pulse  in  flame,  40  Hz,  Rj  10  excitation  (281.537  nm),  50  pm/3  mm  slits. 
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Figure  35.  "OH  LIF”  profile  in  propellant  flame.  AN1  propellant,  0.8  MPa  (100  psig),  0.15 
mJ/pulse  in  flame  (unattenuated),  40  Hz,  10  excitation  (281.537  nm),  50  pm/3  mm  slits. 
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Figure  36.  "OH  LIF"  profiles  in  propellant  flame.  API  propellant,  0.5  MPa  (50  psig),  0.015 
mJ/pulse  in  flame,  40  Hz,  laser  detuned  from  OH  line  (281.48  nm),  50  pm/3  mm  slits. 
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Figure  37.  "OH  LIF"  profile  in  propellant  flame.  HMX1  propellant,  1.5  MPa  (200  psig),  0.12 
mJ/pulse  in  flame,  40  Hz,  laser  detuned  from  OH  line  (281.48  nm),  50  |im/3  mm  slits. 
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Figure  38.  "OH  LIF"  profiles  in  propellant  flame  with  313  nm  filter.  API  propellant,  0.5  MPa 
(50  psig),  0.08  mJ/pulse  in  flame,  40  Hz,  3  pulse  average,  Ri  10  excitation,  30  p.m/3  nmi  slits. 


31 


Distance  above  surface,  mm 


Figure  39.  "OH  LIF"  profiles  in  propellant  flame  with  313  nm  filter.  API  propellant,  0.5  MPa 
(50  psig),  0.08  mJ/pulse  in  flame,  40  Hz,  3  pulse  average,  Rj  5  excitation,  30  p.m/3  mm  slits. 
One  experiment  was  truncated  by  experimental  difficulties  when  the  laser  was  ~2mm  above  the 

surface. 


Distance  above  surface,  mm 


Figure  40.  "OH  LIF"  profiles  in  propellant  flame  with  313  nm  filter.  API  propellant,  0.5  MPa 
(50  psig),  0.08  mJ/pulse  in  flame,  40  Hz,  3  pulse  average,  laser  detuned  from  OH  lines,  30  |im/3 

mm  slits. 
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Distance  above  surface,  mm 


Figure  41.  "OH  LIF"  profiles  in  propellant  flame  with  313  nm  filter.  HMX1  propellant,  1.5  MPa 
(200  psig),  0.08  mJ/pulse  in  flame,  40  Hz,  3  pulse  average,  Rj  5  excitation,  30  p.ni/3  mm  slits. 
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Distance  above  surface,  mm 


Figure  42.  "OH  LIF"  profiles  in  propellant  flame  with  313  nm  filter.  HMX1  propellant,  1.5  MPa 
(200  psig),  0.08  mJ/pulse  in  flame,  40  Hz,  3  pulse  average,  Rj  10  excitation,  30  (im/3  mm  slits. 
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Distance  above  surface,  mm 


Figure  43.  "OH  LIF"  profiles  in  propellant  flame  with  313  nm  filter.  HMX1  propellant,  1.5  MPa 
(200  psig),  0.08  mj/pulse  in  flame,  40  Hz,  3  pulse  average,  laser  detuned  from  OH  lines,  30  p.m/3 

mm  slits. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


This  report  has  described  an  attempt  to  use  OH  LIF  to  measure  temperature  in  high 
pressure  solid  propellant  flames.  The  techniques  used  were  similar  to  those  successfully  employed 
by  Parr  and  Parr  (Ref.  13)  to  measure  temperature  in  atmospheric  pressure  laser-driven  propellant 
"flames."  When  applied  to  higher  pressure  propellant  flames,  however,  the  requirements  of  a  large 
(20  nm)  detection  bandpass  to  obtain  accurate  measurements  (Ref.  28,  50)  created  severe  problems 
with  light  rejection.  The  attenuation  of  the  diagnostic  laser  beam  caused  by  the  optical  thickness  of 
the  propellant  flames  created  very  large  signals  at  the  laser  wavelength  (281  nm),  which  were 
difficult  to  separate  from  the  desired  LIF  signal  (315±10  nm).  Using  the  combination  of  an 
interference  filter  and  a  spectrometer  (300  gr/mm  grating  with  a  3-mm  rear  slit),  wide  bandpass 
detection  with  relatively  high  light  rejection  could  be  accomplished,  although  only  marginally 
successful  results  could  be  obtained  in  an  AP-based  propellant  with  its  high  OH  concentration. 
OH  LIF  temperature  measurements  in  HMX-based  propellant  flames  were  unsuccessful.  Because 
of  the  difficulty  of  the  OH  LIF  temperature  measurements  in  high  pressure  propellant  flames,  other 
temperature  measurement  techniques,  such  as  CARS,  may  be  more  successful  and  should  be 
pursued.  Another  option  is  NO  LIF  temperature  measurements.  NO  has  a  much  different  vibro- 
rotational  structure  from  OH,  so  it  may  be  possible  to  accomplish  the  necessary  wide  bandpass 
detection  (which  would  be  much  less  than  the  20  nm  required  for  OH)  with  the  necessary  light 
rejection. 
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APPENDIX  A.  OH  SPECTROSCOPIC  DATA 
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3311.83 

R2  1 

282.4096 

144 

0.5 

1.5 

126.43 
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OH  1,1  band  Wavelength,  B  coefficient  J"  J'  Rotational  energy, 

ID _ nm  (Ref. 39) _ (Ref.  41) _ cnrr1  (Ref.  39} 


Ql  4 

282.4393 

361 

4.5 

4.5 

355.87 

Q1  4' 

282.4455 

53 

4.5 

3.5 

355.87 

R?  13 

282.5545 

212 

12.5 

13.5 

3344.34 

Pi  2 

282.5802 

225 

2.5 

1.5 

83.70 

R,  14 

282.6058 

211 

14.5 

15.5 

3810.60 

Rl  14’ 

282.6317 

11 

14.5 

14.5 

3810.60 

Ql  5 

282.6674 

377 

5.5 

5.5 

544.82 

Q,  5' 

282.6764 

40 

5.5 

4.5 

544.82 

R2  14 

282.8784 

214 

13.5 

14.5 

2841.51 

Q2  2* 

282.8942 

103 

1.5 

2.5 

187.47 

Q2  2 

282.8984 

294 

1.5 

1.5 

187.47 

Q21* 

282.9225 

286 

0.5 

0.5 

126.12 

a2r* 

282.9225 

144 

0.5 

1.5 

126.12 

Ql  6* 

282.9277 

388 

6.5 

6.5 

769.17 

Q2  3’* 

282.9277 

76 

2.5 

3.5 

288.83 

Q!  6'* 

282.9373 

31 

6.5 

5.5 

769.17 

Q2  3* 

282.9373 

323 

2.5 

2.5 

288.83 

Ri  15 

282.9879 

212 

15.5 

16.5 

4341.70 

Pi  3 

283.0095 

218 

3.5 

2.5 

201.90 

Q2  4' 

283.0257 

58 

3.5 

4.5 

429.23 

Q2  4 

283.0339 

348 

3.5 

3.5 

429.23 

Q2  5' 

283.1748 

45 

4.5 

5.5 

608.16 

Q25* 

283.1838 

366 

4.5 

4.5 

608.16 

p2r 

283.1838 

286 

0.5 

0.5 

126.43 

Ql  7 

283.2222 

397 

7.5 

7.5 

1029.10 

Ql  7* 

283.2351 

25 

7.5 

6.5 

1029.10 

R2  15 

283.2512 

215 

14.5 

15.5 

4371.18 

Q2  6’ 

283.3731 

35 

5.5 

6.5 

824.76 

Q2  6 

283.3843 

379 

5.5 

5.5 

824.76 

P2  2’ 

283.4143 

161 

1.5 

1.5 

187.71 

R,  16* 

283.4172 

213 

16.5 

17.5 

4904.53 

P2  2* 

283.4172 

116 

1.5 

0.5 

187.71 

Pi  4 

283.4630 

215 

4.5 

3.5 

355.09 

Qi  8 

283.5530 

403 

8.5 

8.5 

1324.24 

Q\  8’ 

283.5675 

20 

8.5 

7.5 

1324.24 

Q2  7' 

283.6189 

28 

6.5 

7.5 

1078.47 

Q2  7 

283.6320 

389 

6.5 

6.5 

1078.47 

R2  16* 

283.6736 

215 

15.5 

16.5 

4932.80 

02  2* 

283.6736 

34 

1.5 

0.5 

187.47 

P2  3’ 

283.7094 

115 

2.5 

2.5 

289.01 

P2  3 

283.7134 

150 

2.5 

1.5 

289.01 

R,  17 

283.8984 

213 

17.5 

18.5 

5498.27 

Q2  8' 

283.9108 

23 

7.5 

8.5 

1368.66 

Qt  9* 

263.9234 

407 

9.5 

9.5 

1654.50 

Q2  8* 

283.9234 

397 

7.5 

7.5 

1368.66 

Qi  9' 

283.9378 

16 

9.5 

8.5 

1654.50 

Pi  5 

283.9450 

214 

5.5 

4.5 

543.54 

OH  1 ,1  band 
ID 


P24* 
P24 
R2  17 
02  3 
Q2  9’ 
Q2  9 

Ql  10 

Q:  10' 
18 


P2  5 
Q210* 
Q2  10 
R2  18 
Q1  11 
0,11* 
02  4 
P2  6' 
P2  6 


284.0621 

86 

3.5 

3.5 

284.0681 

169 

3.5 

2.5 

284.1481 

216 

16.5 

17.5 

284.2282 

30 

2.5 

1.5 

284.2471 

19 

8.5 

9.5 

284.2635 

402 

8.5 

8.5 

284.3294 

411 

10.5 

10.5 

284.3470 

14 

10.5 

9.5 

284.4304 

214 

18.5 

19.5 

284.4581 

213 

6.5 

5.5 

284.4670 

66 

4.5 

4.5 

284.4750 

181 

4.5 

3.5 

284.6280 

16 

9.5 

10.5 

284.6461 

407 

9.5 

9.5 

284.6741 

216 

17.5 

18.5 

284.7771 

413 

11.5 

11.5 

284.7967 

12 

11.5 

10.5 

284.8429 

24 

3.5 

2.5 

284.9206 

52 

5.5 

5.5 

284.9304 

188 

5.5 

4.5 

285.0046 

212 

7.5 

6.5 

Rotational  energy, 
cm'1  (Ref.  39) 


429.45 

429.45 

5525.50 
288.83 
1694.85 
1694.85 
2019.53 
2019.53 

6122.51 

767.45 
608.15 
608.15 
2056.46 
2056.46 
6148.33 
2419.05 
2419.05 
429.23 
824.49 
824.49 
1026.69 


Overlapped  by  nearby  lines  and  thus  not  usable.  The  rotational  energy  refers  to  the  lower  state  (v”). 
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APPENDIX  B.  OH  SPECTRAL  FITTING  PROGRAM 

As  a  part  of  this  effort,  a  program  to  create  "synthetic"  OH  spectra  was  constructed.  This 
is  a  limited  effort  compared  to  the  program  created  by  Campbell  (Ref.  11),  who  included  such 
effects  as  self-absorption  and  a  fitting  routine  to  simulate  OH  emission  spectra  from  high  pressure 
solid  propellant  flames. 

The  program  calculates  the  line  locations  and  intensities  of  the  12  branches  of  the  OH  A2!  - 
X2I1  system.  The  branches  are  listed  in  Table  2.  The  line  locations  were  originally  calculated 
using  the  standard  spectroscopic  equations  for  the  OH  transition  terms,  e.g., 

Fj(v,N)=Bv(v)*N*(N+l)-Dv(v)*N2*(N+l)2+R*(N+0.5)+vibrational  terms. 

However,  this  yielded  unsatisfactory  results  for  the  line  locations,  so  tabular  values  for  the  terms 
were  used  in  the  program  (Ref.  39,48).  Then  the  wavenumber  of  the  line  was  calculated  as  the 
differences  of  the  terms,  as  shown  in  the  table.  The  relative  intensity  of  a  given  line  was  calculated 
as  E(v,N)=A*g*BF*BG;  where  A  in  the  Einstein  coeffcient  for  stimulated  emission;  g  is  the 
degeneracy  (2J+1);  BF  is  the  rotational  Boltzmann  factor  (BF=exp(-hcF/kTf));  and  BG  is  the 
vibrational  Boltzmann  factor  (BG=exp(-hcG/kTv)).  The  Einstein  coefficients  and  vibrational  terms 
were  taken  from  tables  (Ref.  41,39).  The  vibrational  and  rotational  temperatures  (Tv  and  Tr, 
respectively)  were  inputs  to  the  program.  When  the  lines  and  intensities  were  calculated,  the 
program  "scanned"  a  triangular  slit  function  over  the  lines  to  approximate  the  experimental  data 
taken  with  spectrometers.  For  this  initial  program,  the  width  of  the  lines  was  assumed  to  be  very 
small  compared  to  the  slit  function,  with  thq  lines  treated  as  delta  functions  at  the  line  center.  The 
width  (FWHM)  of  the  slit  function,  the  number  of  steps  in  the  "scan,"  and  the  incremental  distance 
between  steps  were  inputs  to  the  program.  To  model  emission  spectra  or  some  aspects  of  the  LIF 
spectra,  the  values  of  A,  g,  BF,  and  BG  were  calculated  for  the  excited  state  (A2E).  For  example, 
a  synthetic  OH  emission  spectrum  with  Tv=Tr=2500  K  is  shown  in  Figure  44.  The  scanning 
parameters  were  selected  to  simulate  emission  spectra  taken  with  a  Reticon  diode  array  detector 
(e.g..  Figure  18).  The  slit  function  has  a  FWHM  width  of  0.125  nm,  1500  data  points,  with 
0.0175  nm  between  points  starting  at  305  nm.  The  match  between  experimental  and  synthetic 
spectra  is  fairly  good.  To  simulate  a  low  pressure  LIF  spectrum  where  there  is  little  collisional 
redistribution  before  radiation,  the  program  could  be  constrained  to  examine  only  lines  directly 
coupled  to  the  laser  line.  An  example  of  such  a  calculation  is  shown  in  Figure  45.  It  is  obvious  by 
comparing  Figures  45  and  2  that  a  great  deal  of  rotational  and  vibrational  redistribution  is 
occurring.  If  one  assumes  total  rotational  thermalization  and  limited  vibrational  thermalization,  the 
OH  model  spectrum  shown  in  Figure  46  is  obtained.  It  is  apparent  from  this  that  the  LIF  spectrum 
of  Figure  2  is  partly  thermalized  by  collisions.  Detailed  models  of  flame  LIF  experiments  can  be 
used  to  calculate  a  temperature  from  a  spectrum  similar  to  Figure  2  (Ref.  51,52),  although  the 
models  are  complex  and  the  detailed  rate  data  needed  is  not  totally  reliable. 


45 


Table.  OH  model  data 


Transition _ 1  N' 


loFiTTloriiiHifi 


N"-l 

P2=F2(N’)-f2(N")  N"-l 

Qi-FjOO-fj'CN")  N" 

Q2*F2(N,)-f2,(N")  N" 

R1=F1(N')-f1(N")  N"+l 

R2=F2(N’)-f2(N")  N"+l 

Satellite  transitions 


012=02=F1(N')-f2,(N'')  N"-2 

Pl2=P2’=Fi(N,)-f2(N")  N"-l 

Ql2=Q2’=Fl(N’H2'(N")  N" 

Q21=Ql’=F2(N,)-fr(N")  N" 

R21=R1'=F2(N’)-f1(N")  N"+l 

S2i=Si=F2(N')-f1'(N")  N"+2 


j" _ y _ 

N"+0.5  J"-l=N"-0.5=N'+0.5 

N"-0.5  J"-l=N"-1.5=N'-0.5 

N"+0.5  J"=N"+0.5=N’+0.5 

N"-0.5  J"=N"-0.5=N'-0.5 

N"+0.5  J"+l=N"+1.5=N'+0.5 

N"*0.5  J"+l=N"+0.5=N'-0.5 


N”-0.5  J"-l=N"-1.5=N,+0.5 

N"-0.5  J"=N"-0.5=N’+0.5 

N"-0.5  J"+l=N"+0.5=N'+0.5 

N"+0.5  J"-l=N"-0.5=N'-0.5 

N"+0.5  j"=N"+0.5=N’-0.5 

N"+0.5  J"+l=N"+1.5=N'-0.5 


OH  emission  intensity 


50 
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Figure  46.  OH  model  spectrum.  Lines  directly  coupled  to  OH  1,0  R!  5. 
Lines  are  1,1  5,  Qi  6,  Q12  6  (weak),  Pj  7,  and  P22  7  (very  weak). 


Figure  47. 


t.  Same  as  Figure  44  but  0,0  A  coefficients  divided  by  50. 


Program  listing:  (run  on  Macintosh  computer  using  Absoft  Fortran  Compiler  Version  2.3) 


PROGRAM  BOLTZM 
C 

C  CALCULATES  SYNTHETIC  SPECTRA  IN  MOLECULE 
C  AS  A  FUNCT!  ON  OF  TEMPERATURE .  CURRENTLY  ONLY  FOR 
C  0,0  AND  1,1  BANDS.  BASED  ON  BOLTZMANN  DISTRIBUTION. 

C  ALLOW  FOR  V=l-4  AND  N=l-40. 

C 

DIMENSION  VZERO ( 4 ,  4 ) , G1 ( 5 ) , G2 ( 5 ) 

DIMENSION  ENET (2000) , EINCR (2000) , WSCAN (2000) , ESUM(2000) 
DIMENSION  Fll (4,  40)  ,F12 (4, 40) , F21 (4, 40) , F22 (4, 40) 

DIMENSION  VP1 (4, 40) ,  VP2 (4, 40) , VQ1 (4, 40) , VQ2 (4, 40) 

DIMENSION  VQ12 (4, 40) ,VQ21 (4, 40) , VR21 (4, 40) ,VP12 (4, 40) 

DIMENSION  VR1 (4, 40) ,VR2 (4,  40) 

DIMENSION  F21P (4, 40) , B1 (4) , F22P (4,  40) 

DIMENSION  BF1 (4, 40) , BF2 (4, 40) , BG  (4) 

DIMENSION  ENR1 (4,  40) ,  ENR2 (4, 40)  ,ENQ1 (4,  40)  ,ENQ2 (4,40) 
DIMENSION  ENP1 (4, 40) ,ENP2 (4,40), ENR21 (4,40) , ENP12 (4,40) 
DIMENSION  ENQ12 (4,40), ENQ21 (4,  40) 

DIMENSION  WAVQ1 (4,40), WAVQ2 (4,40), WAVP1 (4,40) 

DIMENSION  WAVP2 (4, 40) , WAVR1 (4, 40) , WAVR2 (4,  40) 

DIMENSION  WAVR21 (4, 40) ,WAVP12 (4, 40) ,WAVQ12 (4, 40) ,WAVQ21 (4, 40) 
DIMENSION  ASPP1 (4, 40) , ASPP2 (4, 40) ,ASPQ1 (4, 40) 

DIMENSION  ASPQ2 (4, 40) , ASPR1 (4, 40) , ASPR2 (4,40) 

DIMENSION  ASPQ12 (4,40) , ASPQ21 (4,40) , ASPR21 (4, 40) , ASPP12 (4,40) 
DIMENSION  GR1 (25) ,GR2 (25) ,GQ1 (25) ,GQ2 (25) ,GP1 (25) ,GP2 (25) 
DIMENSION  GP12 (25) ,GR21 (25) , GQ12 (25) ,GQ21 (25) 

DIMENSION  G012 (25) , GS21 (25) 

INTEGER* 2  N, COLUMNS 
INTEGER* 2  NDATA, DELTV 
C 

OPEN (UNIT=23,  FILE='OH  START',  FORM= ' FORMATTED *, STATUS= ' OLD  * ) 
READ  (23,250)  TROT, TVIB,WINIT, RES, WINCR, NDATA, DELTV 
250  FORMAT (3F9. 3,F6. 3,F6. 4,15,14) 

C 

C  CAN  WRITE  MOLECULE  INFO  IF  DESIRED 
C 

WRITE (9, 300) 

300  FORMAT  ('MOLECULE  IS  OH') 

WRITE (9, 400)  TROT, TVIB 

400  FORMAT  (' BOLTZMANN  DISTRIBUTION  FOR  TROT, TVIB= ', F9 . 3, F9 . 3) 
WRITE (9,450)  WINIT, RES, WINCR, NDATA,  DELTV 
450  FORMAT ('WINIT=',F9. 3, 5X, 'RES=' ,F6. 3, 5X, 'WINCR=' , F6. 4, /, 

1  ' NDATA= ' , I 4 , 5X, ' DELTV= ’,14) 

CLOSE (UNIT=23) 

OPEN (UNIT=23,  FILE='OH  LINES',  FORM= ' FORMATTED ' , STATUS= ' NEW ' ) 

C 

C 

C  NOTE ! ! ! ! !  NEED  DUMMY  VARIABLES  SlfoCE  ARRAYS  CANT  HAVE  ZEROS ! 

C  SO  NV1=NV2=1  FOR  0,0  BAND,  ETC,  AND  Nl-0,  SO  USE  N1  IN  CALCS 

C 

C  NEED  DATA  FOR  EINSTEIN  A  COEFFICIENTS  FOR  SPONTANEOUS 

C  EMISSION  (ASP** (NV1,N)=A  COEFF  FOR  BRANCH  **) 

C  FROM  CHIDSEY  AND  CROSLEY  (JQSRT  23:187-199  (1980)) 

C 

ASPP1 (1,1) =1000 


48 


onoooo  oooocj  co  oo  ooo  on  non  non  ooooo  o  o  on 


ASPP1 (1,2) =671 


(data  omitted  to  save  paper) 

ASPR2  (2, 27) =60 
ASPR2 (2, 28) =53 

DATA  FROM  HUBER  AND  HERZBERG 
T=ELECTRONIC  TERM  FOR  OH  A  STATE 

Tl=32684 . 1 

DATA  FROM  DIEKE  AND  CROSSWHITE 
FOR  G (V) “VIBRATIONAL  (ROTATIONLESS)  TERM  ENERGY 
G1=G'  (A  STATE),  G2=G"  (GROUND  STATE) 

G1 (1) =1567. 68 
G1 (2)=G1 (1) +2988. 60 
G1  (3)=G1 (2) +2792. 92 
G1 ( 4 ) =G1 (3) +2593.36 
G2 (1)=1847.052 
G2 (2) =G2 (1) +3569.59 
G2  (3)=G2 (2) +3403. 97 


DO  2600  NV1=1, 2 

HERE  IS  WHERE  THE  BAND  SEQUENCES  ARE  INPUT— 0, 0->NVl=NV2- 
THUS  DELTV=0, +1,-1  FOR  0>0; 0,1;  AND  THE  1,0  SEQUENCES 
NV2=NV1-DELTV 


VZERO (NV1 , NV2 ) =T1 +G1 (NV1 ) -G2 (NV2 ) 

N1,N2  ARE  ACTUAL  N 

N1=NV1-1 

N2=NV2-1 


WRITE (9, 350)  N1,N2 

50  FORMAT  (/, 'CALCULATING' , 3X, 13, 13, 3X, 'LINES' ) 

WRITE (23, 375)  N1,N2 
75  FORMAT  (/, 'BAND  IS' , 13, ' , ' , 13) 

VIBRATIONAL  BOLTZMANN  FACTOR=BG (NV1 ) 
hc/k=l. 43836  [K/WAVENUMBER] 

BG (NV1 ) =EXP (-1 . 43836*G1 (NV1) /TVIB) 

DATA  INPUT  FORMAT  F2=f,  F1=F 

F21=fl  COMPONENT,  F21P=fl',  ETC.,  FUNCTION  OF  N  AND  V 
F (N, NV)  LISTED  BELOW.  F  AND  f  IN  TABLES  INCLUDES 
VIBRATIONAL  AND  ELECTRONIC  ENERGY 

Fll  (1,1)=32474.62 
Fll  (1, 2) =32542 . 56 
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(data  omitted  to  save  paper) 

F22P (2, 27) =16069. 01 
F22P (2, 28) =16915. 75 
C 

DO  2500  N=l, 25 
C 

C  N=N",  NP=N '  FOR  P  BRANCH,  ETC. 

C 

NP=N-1 

NR=N+1 

NQ=N 

NO=N-2 

NS=N+2 

C 

C  G=DEGENERACY=2J'+1  IN  EMISSION 
C 

GR1 (N) =2* (N+l . 5) +1 
GR2 (N) =2* (N+0. 5) +1 
GQ1 (N) =2* (N+0 . 5) +1 
GQ2 (N) =2* (N-0 . 5) +1 
GP1(N) =2* (N-0 . 5) +1 
GP2 (N) =2* (N-l . 5) +1 
G012 (N) =2* (N-2 . 5) +1 
GP12 (N) =2* (N-l . 5) +1 
GQl2(N)=2*(N-0.5)+l 
GQ21  (N)  =2*  (N+0 . 5)  +1 
GR21  (N)  *2*  (N+1.5J+1 
GS21  (N)  *2*  (N+2-.  5)  +1 
C 

IF (NP.EQ. 0)  GO  TO  650 
C 

VP1 (NV1, N) =F11 (NV1, NP) -F21 (NV2,  N) 

VP 2 (NV1,N)=F12 (NV1, NP) -F22 (NV2,N) 

VP12 (NV1,N)=F11 (NV1,NP) -F22 (NV2,N) 

C 

C  INPUT  DATA  FOR  0,0  AND  1,1  BANDS  FOR  PI  AND  P2,  N=1 
C 

650  VP1 (1,1) =32440. 60 
VP2 (1, 1) =0 
VP1 (2,1) =31860. 78 
VP2 (2, 1)=0 
VP12 (1,1) =32314. 19 
VP12 (2,1)=31733.71 
C 

C  IGNORE  O  AND  S  BRANCHES  FOR  NOW 
C 

VQ1 (NV1,N) =F11 (NV1,NQ) -F21P (NV2,N) 

VQ2 (NV1 , N) =F12 (NV1 , NQ) -F22P (NV2 , N) 

VR1 (NV1,N)»F11 (NV1, NR) -F21 (NV2,N) 

VR2 (NV1, N) =F12 (NV1,NR) -F22 (NV2,N) 

VQ12 (NV1,N) =F11 (NV1, NQ) -F22P (NV2,N) 

VQ21 (NV1, N) =F12 (NV1,NQ) -F21P (NV2,N) 

VR21 (NV1 , N) =F12 (NV1 , NR) -F21 (NV2 , N) 

C 

C  BF=ROTATIONAL  BOLTZMANN  FACTOR 

C  SUBTRACT  VIBRATIONAL  ENERGY  TO  OBTAIN  ROTATION  ONLY 
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c 

BF1  (NV1,  N)  =EXP  (-1*1 . 43836*  (Fll  (NV1,  N)  -Tl-Gl  (NV1)  )  /TROT) 
BF2 (NV1,N)=EXP (-1*1. 43836* (F12 (NV1,N) -Tl-Gl (NV1) ) /TROT) 
C 

ENR1  (NV1,  N)  =ASPR1  (NV1,  N)  *BF1  (NV1,  N)  *BG  (NV1)  *GR1  (N) 

ENR2  (NV1 ,  N)  =ASPR2  (NV1,  N)  *BF2  (NV1 ,  N)  *BG  (NV1 )  *GR2  (N) 

WAVR1 (NV1 , N) =1 . 0E07/ (VR1 (NV1, N) *1 . 00029) 

WAVR2 (NV1, N) =1 . 0E07/ (VR2 (NV1,N) *1.00029) 

C 

ENQ1 (NV1 , N) =ASPQ1 (NV1 , N) *BF1 (NV1 ,N) *BG (NV1 ) *GQ1 (N) 

ENQ2 (NV1,N) =ASPQ2 (NV1,N) *BF2 (NV1,N) *BG (NV1) *GQ2 (N) 
WAVQ1 (NV1,N)=1.0E07/ (VQ1 (NV1,N) *1.00029) 

WAVQ2 (NV1 , N) =1 . 0E07 / (VQ2 (NV1 ,N) *1.00029) 

C 

ENP12 (NV1 , N ) =ASPP1 2 (NV1,N) *BF1 (NV1 , N) *BG (NV1 ) *GP12 (N)  . 
ENR21 (NV1 , N) =ASPR21 (NV1 , N) *BF2 (NV1 ,N) *BG (NV1 ) *GR21 (N) 
WAVP12 (NV1 , N) =1 . 0E07/ (VP12 (NV1 , N) *1 . 0002 9 ) 

WAVR21 (NV1,N)=1.0E07/(VR21 (NV1, N) *1 . 00029) 

C 

ENQ12  (NV1 ,  N)  =ASPQ12  (NV1 ,  N)  *BF1  (NV1 ,  N)  *BG  (NV1 )  *GQ12  (N) 
ENQ21  (NV1,  N)  =ASPQ21  (NV1,N)  *BF2  (NV1,N)*BG  (NV1)  *GQ21  (N) 
WAVQ12 (NV1, N) =1 . 0E07/ (VQ12 (NV1, N) *1 . 00029) 

WAVQ21 (NV1, N) =1 . 0E07/ (VQ21 (NV1,N) *1.00029) 

C 

C 

ENP1 (NV1 , N) =ASPP1 (NV1,N) *BF1 (NV1,N) *BG (NV1) *GP1 (N) 

ENP2 (NV1, N) =ASPP2 (NV1, N) *BF2 (NV1, N) *BG (NV1) *GP2 (N) 

C 

IF (VP1 (NV1,N) .EQ.0.0)  GO  TO  890 
IF (VP2 (NV1 , N) . EQ .0.0)  GO  TO  890 

C 

WAVP1 (NV1 , N) =1 . 0E07 / (VP1 (NV1,N) *1.00029) 

WAVP2 (NV1, N) =1 . 0E07/ (VP2 (NV1,  N) *1 . 00029) 

C 

890  CONTINUE  • 

C 

WRITE (23, 895)  N, VR1 (NV1 , N) , ENR1 (NV1, N) , VR2 (NV1, N) , 

1  ENR2 (NV1,N) 

895  FORMAT  ( 'N=' ,  12, IX,  'VR1==' ,  F10 . 4, 2X,  ' ENR1= ' ,  F10 . 4 , 2X, 

1  'VR2=',F10.4,2X, 'ENR2=',F10.4) 

WRITE (23,896)  N,VQ1 (NV1,N) , ENQ1 (NV1,N) ,VQ2 (NV1,N) , 

1  ENQ2 (NV1,N) 

896  FORMAT  ('N=',  12,  IX,  ’  VQ1=*  ,F10 . 4,  2X,  'ENQ^'^IO^^X, 

1  •VQ2=,,F10.4,2X/ 'ENQ2=',F10.4) 

WRITE (23, 897)  N,VP1 (NV1,N) ,ENP1 (NV1,N) ,VP2 (NV1,N) , 

1  ENP2 (NV1,N) 

897  FORMAT ('N=’, 12, IX, *VP1=' ,F10 . 4, 2X, 'ENPl=' ,F10. 4, 2X, 

1  'VP2=',F10.4,2X, ,ENP2=',F10.4) 

WRITE (23, 898)  N,VR21 (NV1,N) ,ENR2l (NV1,N) ,VP12 (NV1,N) , 

1  ENP12 (NV1, N) 

898  FORMAT ('N=', 12, IX, ' VR21=' , F10 . 4, 2X, ' ENR21=' ,F10. 4, 2X, 

1  'VP12=',F10.4,2X, 'ENP12=',F10.4) 

WRITE (23, 899)  N, VQ12 (NV1,N) , ENQ12 (NV1,N) , VQ21 (NV1,N) , 

1  ENQ21 (NV1,N) 

899  FORMAT ('N=', 12, IX, ' VQ12=' , F10 . 4, 2X, 'ENQ12=' ,F10 . 4,  2X, 

1  ' VQ21=' , F10 . 4, 2X, 'ENQ21=' ,F10 . 4) 

C 

C 
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2500  CONTINUE 
2600  CONTINUE 

WRITE  (9,2700) 

2700  FORMAT  (/, 'LINE  CALCULATIONS  COMPLETE  -  SCANNING  UNDERWAY',/) 
CLOSE  (UNIT=23) 

OPEN (UNIT=23,  FILE= • DATA • ,  FORM= 1 FORMATTED ' , STATUS= ’ NEW 1 ) 

C 

C  NOW  START  SCANNING  OVER  LINES 
C 

DO  3000  1=1 , NDATA 
3000  WSCAN(I) =0.000 
WSCAN ( 0 ) =WINIT 
DO  8000  L=l, NDATA 
DO  3100  12=1, NDATA 
3100  EINCR (12) =0.000 

DO  3200  13=1, NDATA 
3200.  ENET  (13)  =0.000 
M=L-1 

WSCAN  (L )  =WSCAN  (M)  +WINCR 

C 

C  NOW  FOR  R  LOOPS 
C 

DO  6700  NV1=1 , 2 
C 

DO  4000  N=l, 25 
WDELT=WSCAN (L) -WAVR1 (NV1 ,  N) 

DELTA-ABS (WDELT) 

IF (DELTA. GT. RES)  GO  TO  3900 
FRAC= (RES-DELTA) /RES 
GO  TO  3950 
3900  FRAC=0 .000 
3950  EINCR (N) =FRAC*ENR1 (NV1,  N) 

ENET (L) =ENET (L) +EINCR (N) 

4000  CONTINUE 

DO  4100  N=l, 25 
WDELT=WSCAN (L) -WAVR2 (NV1 ,  N) 

DELTA=ABS (WDELT) 

IF (DELTA. GT. RES)  GO  TO  4025 
FRAC= (RES-DELTA) /RES 
GO  TO  4050 
4025  FRAC=0 .000 
4050  EINCR (N) =FRAC*ENR2 (NV1 , N) 

ENET (L) =ENET (L) +EINCR (N) 

4100  CONTINUE 
C 

C  NOW  FOR  P  LOOPS 
C 

DO  4250  N=l, 25 
WDELT=WSCAN (L) -WAVP1 (NV1 , N) 

DELTA=ABS (WDELT) 

IF (DELTA. GT. RES)  GO  TO  4220 
FRAC= (RES-DELTA) /RES 
GO  TO  4225 
4220  FRAC=0 .000 
4225  EINCR (N) =FRAC*ENP1 (NV1,N) 

ENET (L) =ENET (L) +EINCR (N) 

4250  CONTINUE 

DO  4300  N=l, 25 
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WDELT=WSCAN (L) -WAVP2 (NV1, N) 
DELTA=ABS (WDELT) 

IF (DELTA. GT. RES)  GO  TO  4275 
'  FRAC= (RES-DELTA) /RES 

GO  TO  4280 
4275  FRAC=0.000 

"  4280  EINCR(N)=FRAC*ENP2(NV1,N) 

ENET (L) =ENET (L) +EINCR (N) 

4300  CONTINUE 
C 

C  NOW  FOR  Q  LOOPS 
C 

DO  4500  N=l, 25 
WDELT=WSCAN (L) -WAVQ1 (NV1,N) 
DELTA=ABS (WDELT) 

IF (DELTA. GT. RES)  GO  TO  4425 
FRAC= (RES-DELTA) /RES 
GO  TO  4450 
4425  FRAC=0 .000 
4450  EINCR  (N) =FRAC*ENQ1 (NV1 , N) 
ENET (L) =ENET (L) +EINCR (N) 

4500  CONTINUE 

DO  4600  N=l, 25 
WDELT=WSCAN (L) -WAVQ2 (NV1 , N) 
DELTA=ABS (WDELT) 

IF (DELTA. GT. RES)  GO  TO  4525 
FRAC= (RES-DELTA) /RES 
GO  TO  4550 
4525  FRA C-0.000 
4550  EINCR (N) -FRAC*ENQ2 (NV1 , N) 
ENET (L) =ENET (L) +EINCR (N) 

4600  CONTINUE 
C 

C  NOW  FOR  Q21,Q12  LOOPS 
C 

DO  5500  N=l, 25 
WDELT=WSCAN(L) -WAVQ21 (NV1,N) 
DELTA-ABS (WDELT) 

IF (DELTA. GT. RES)  GO  TO  5425 
FRAC= (RES-DELTA) /RES 
GO  TO  5450 
5425  FRAC=0 .000 
5450  EINCR (N) =FRAC*ENQ21 (NV1,  N) 
ENET (L) =ENET (L) +EINCR (N) 

5500  CONTINUE 

DO  5600  N»l,25 
WDELT=WSCAN (L) -WAVQ12 (NV1 , N) 
DELTA=ABS (WDELT) 

IF (DELTA. GT. RES)  GO  TO  5525 
FRAC- (RES-DELTA) /RES 
GO  TO  5550 
5525  FRAC»0.000 

5550  EINCR (N) »FRAC*ENQ12 (NV1 , N) 
ENET (L) =ENET (L) +EINCR (N) 

5600  CONTINUE 
'  C 

C  NOW  FOR  R21, P12  LOOPS 
C 
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6425 

6450 

6500 


6525 

6550 

6600 

C 

6700 

C 

C 

C 

7895 

C 

C 


7900 

7950 

8000 

8100 

C 

C 

C 

8200 

8300 

C 

C 

C 

C 

C 

c 

C 


DO  6500  N=l, 25 
WDELT=WSCAN (L) -WAVR21 (NV1,N) 

DELTA=ABS (WDELT) 

IF (DELTA. GT. RES)  GO  TO  6425 
FRAC= (RES-DELTA) /RES 
GO  TO  6450 
FRAC=0 .000 

EINCR (N) =FRAC*ENR21 (NV1,  N) 

ENET (L) =ENET (L) +EINCR (N) 

CONTINUE 

DO  6600  N=l, 25 

WDELT^WSCAN (L) -WAVP12 (NV1 , N) 

DELTA=ABS (WDELT) 

IF (DELTA. GT. RES)  GO  TO  6525 
FRAC= (RES-DELTA) /RES 
GO  TO  6550 
FRAC=0.000 

EINCR  (N)  “FRAC * ENP 1 2  (NV1,  N) 

ENET (L) =ENET (L) + EINCR (N) 

CONTINUE 

CONTINUE 

NOW  WRITE  DATA  TO  FILE 

WRITE (23, 7895)  WSCAN(L) , CHAR (9) ,ENET(L) 
FORMAT  (F8.3,A1,F8.2) 


IF  (L.EQ.100)  GO  TO  7900 

IF  (L.EQ.200)  GO  TO  7900 

IF  (L.EQ.300)  GO  TO  7900 

IF  (L.EQ.400)  GO  TO  7900 

IF  (L.EQ.500)  GO  TO  7900 

IF  (L.EQ.600)  GO  TO  7900 

IF  (L.EQ.750)  GO  TO  7900 

IF  (L.EQ.1000)  GO  TO  7900 

IF  (L.EQ.1250)  GO  TO  7900 

IF  (L.EQ.1500)  GO  TO  7900 

GO  TO  8000 

WRITE  (9,7950)  L,NDATA 

FORMAT  (14, '/', 14, 2X, 'DATA  POINTS') 

CONTINUE 
WRITE  (9,8100) 

FORMAT  (/, 'SCANNING  COMPLETE  -  POST-PROCESSING  UNDERWAY',/) 


CALL  FILECR ( ' DATA ' ) 

CLOSE  (UNIT=23) 

END 


This  subroutine  changes  the  File  whose  name  is  passed 
to  the  type  and  creator  of  a  KALEIDAGRAPH  text 
file.  Created  and  modified  by  F.  Frederick  II  and 
J.  Amfahr. 
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ooo  ooo  ooo  non  ooo 


SUBROUTINE  FILECR(  NAME  ) 

C 

character* (*)  NAME 
INTEGER* 4  toolbx 
CHARACTER* 256  FILENAME 
INTEGER  PTR 

PARAMETER  (PTR=Z ' COOOOOOO ' ) 

C 

include  XP  20: FORTRAN: INCLUDE  FILES: file. inc  . 

include  XP  20: FORTRAN: INCLUDE  FILES : params. inc 

change  the  file  name  to  a  PASCAL  string 

FILENAME  «  char (len (TRIM (  NAME) )) //TRIM (NAME  ) 

set  the  toolbox  parameters  to  call  PBGETFILEINFO 

ionameptr  =  toolbx (PTR,  FILENAME  ) 
iocampletion  =  0 
iofdirindex  =  -1 
iocampletion  =  0 
iovrefnum  =  0 
iofversnum  =  0 

toolbox  routine  to  get  file  info 

CALL  toolbx (PBGETFILEINFO,  toolbx (PTR,  params) ) 

Kaleidagraph  text  file  type  and  creator 

fdtype  =  'TEXT' 
fdcreator  =  'QKPT' 

toolbox  routine  to  set  file  info 

CALL  toolbx (PBSETFILEINFO,  tOOlbx(PTR,  params)) 

RETURN 
END 
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